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f~^ , ABSTRACT. A description of Lagrangian and Hamiltonian formalisms naturally 

^D I arisen from the invariance structure of given nonlinear dynamical systems on the infinite- 

dimensional functional manifold is presented. The basic ideas used to formulate the 
canonical symplectic structure are borrowed from the Cartan's theory of differential sys- 
tems on associated jet-manifolds. The symmetry structure reduced on the invariant 
submanifolds of critical points of some nonlocal Euler-Lagrange functional is described 
thoroughly for both differential and differential discrete dynamical systems. The Hamil- 
tonian representation for a hierarchy of Lax type equations on a dual space to the Lie 
algebra of integral-differential operators with matrix coefficients, extended by evolutions 
for eigenfunctions and adjoint eigenfunctions of the corresponding spectral problems, is 
obtained via some special Backlund transformation. The connection of this hierarchy 
with integrable by Lax spatially two-dimensional systems is studied. 
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One of the fundamental problems in modern theory of infinite-dimensional dynam- 



(N, 

^ : 1 Introduction 

^, 

O . ical systems is that of an invariant reduction them upon some invariant subman- 

O 



ifolds with enough rich mathematical structures to treat their properties analyti- 
cally. The first approaches to these problems were suggested still at the late times 
of the preceding century, in the classical oeuvres by S.Lie, J.Liouville, J.Lagrange, 
V.R.Hamilton, J.Poisson and E.Cartan. They introduced at first the important 
concepts of symmetry, conservation law, symplectic, Poisson and Hamiltonian 
structures as well invariant reduction procedure, which appeared to be extremely 
/S ' useful for proceeding studies. These,, notions were widely generalized further by 

5_] ■ „ Ippft Wein,ffars Lax -^ ^ -^ 

Marsden and Weinstein 1 3,4, 20,TrLa?, IL3 „ Bogoyavlensky and Novikov 



louriau 



l[VJ, as well by many other researchers l[8, 10, IL hZ, 131. It ^ems worthwhile to 
mention here also the recent enough studies in |:^:^,':^5, 'zA, 5y, 26, '2{, 'zh, '2^], where 
the special reduction methods were proposed for the integrable nonlinear dynam- 
ical systems on both functional and operator manifolds. In the present paper we 
describe in detail the reduction procedure for infinite dimensional dynamical sys- 
tems upon the invariant set of critical points of some global invariant functional. 



The method uses the Cartan's differential-geometric treating of differential ideals 
in Grassmann algebra over the associated jet-manifold. As one of main results, we 
show also that both the reduced dynamical systems and their symmetries, generate 
the Hamiltonian flows on the invariant critical submanifolds of local and nonlocal 
functionals with respect to the canonical symplectic structure upon it. These re- 
sults are generalized for the case of differential-difference dynamical systems being 
given on discrete infinite-dimensional manifolds. The direct procedure to construct 
the invariant Lagrangian functionals for a given apriori Lax-type integrable dynam- 
ical system is presented for both the differential and the differential-difference cases 
of the manifold M. Some remarks on the Lagrangian and Hamiltonian formalisms, 
concerned to infinite-dimensional dynamical systems with symmetries are given. 
The Hamiltonian representation for a hierarchy of Lax type equations on a dual space 
to the Lie algebra of integral-differential operators with matrix coefficients, extended by 
evolutions for eigenfunctions and adjoint eigenfunctions of the corresponding spectral 
problems, is obtained via some special Backlund transformation. The connection of this 
hierarchy with integrable by Lax spatially two-dimensional systems is studied. 

2 General setting 

We are interested in treating a given nonlinear dynamical system 

du/dt = K[u], (2.1) 

with respect to an evolution parameter t G M on an infinite-dimensional functional 
manifold M C C''°°^(]R; M™), possessing two additional ingredients: a homoge- 
nous and autonomous conservation law C G D{M) and a number of homogenous 
autonomous symmetries du/dtj = Kj[u], j = l,k, with evolution parameters 
tj G M. The dynamical system (2.1) is not-supposed to be in general Hamilto- 
nian, all the maps K,Kj : M — > T{M), j = l,k, being considered smooth and 
well-defined on M. 

To pose the problem to be discussed further more definitely, let us involve the 
jet-manifold J*^°°) (M; M™) locally isomorphic to the functional manifold M men- 
tioned above. This means the following: the vector field (2.1) on M is cornpletelv 
equivalent to that on the jet-manifold J*^°°)(]R; M™") via the representation l[l, z] 

{M3u^ K[u]) ^ ( {K{u, u^'\ ..., m("+^)) ^ {x; u, u^'\ ..., u^"^^) G J^°^\R; M™)), 

(2.2) 
where n G Z+ is fixed, x G M is the function parameter of the jet-bundle J(-) (M; M™) - 
M, and vr is the usual projection on the base M. Let us allow also that the smooth 
functional C G D{M) is a conservation law of the dynamical system (2.1), that is 



dC/dt = along orbits of (2.1) for all t G M. Due to the jet-representation (2.2) 
we can write the density of the functional C G D{M) in the following form: 

C= I dxC[u], (2.3) 

with MxR'" 9 [x]u] ^ (a;;M,MW,...,M(^+i)) G j(^+^)(K; K™) being the stan- 
dard jet-mapping and a number N G Z_|_ fixed. Besides, the functional (2.3) will be 
assumed to be non-degenerate in the sense that Hessian of C : J*^^+^^(R; R™)— i> R 

has nonvanishing determinant: det\\ — QuiN+i)Qu(r^+i) II ¥" 0- 

3 Lagrangian reduction 

Consider now the set of critical points M„ C M of the functional C G D{M): 

Mm = {u e M : grad C[u] = 0} (3.1) 

5C{u, . . . ,u^^^^^)/5u - the Euler variational 
7 the manifold Mjy C M is smoothly imbedded 
well-defined one due to the condition HessC ^ 0. Besides, the manifold Mjv 
appears to be invariant with respect to the dynamical system (2.1). This means 
in particular that the Lie-derivative of any vector field X : M —^ T{M), tangent 
to the manifold Mn with respect to the vector field (2.1), is again tangent to Mm, 
that isthe foUowung implication 

X[u] G T^{Mm) => [K,X][u] G T^{Mm) (3.2) 

holds for all u G M^. Here we are in a position to begin with a study of the intrinsic 
structure of the manifold Mjv C M. withinthereometric Cartan's theory developed 
on the jet-manifold j(~)(R; R'") pT^T^'Sptet us define an ideal /(^) C A(J(°°)), 
generated by the vector one-forms C = du^^^ — u'^^'^^^dx, j G Z_|_, which are 
canceled the vector field d/dx on the jet-manifold J*^°°)(R; R™) : 

^^d'^ = 0, je Z+, (3.3) 



where, due to (2.2), grad C[u 
derivative. As proved by Lax 



dx 



where a; G R belongs to the jet-bundle base, i_d_- the intrinsic derivative along 

dx 

the vector field 

dx dx ^ ' du^^^ 

iez+ 

where (.,.) is the standard scalar product in R™. The vector field (2.1) on the 
jet-manifold J*^°°^ (R; R) has the analogous representation: 

|-|^5:<*-.al^>. (3.4) 



where, by definition, K^^' := -^K, j = Z+. The problem arises: how to build the 
intrinsic variables on the manifold Mn C M from the jet-manifold coordinates on 
J(°°)(M;M™))? 

To proceed to the solution of the problem above, let us study the 1-form dC = 
^i^j{oo)^]]j. ^y\j g^g Qj^g defined on the submanifold M^ C M. We have the following 
chain of identities in the Grassmann subalgebra A(J'^^^+^^(M; M™)) : 

N 

dC = d{ij_Cdx) = dij_{Cdx+y (pj,]R)) (3.5) 

dx dx ' ^ 

j=0 

N N 

= {dijL + i_d_d) ■ {Cdx + y^{Pj, ^j )) — ijLd{Cdx + /^{PjiCj )); 

i=o " j=o 

where Pj : J*^^^+^)(]R; M™") -^ W^,j = 0,N, are some until not definite vector- 
functions. Requiring now that 2-form d{Cdx + Ylj=o{Pj^^j )) "^^ not depend on 
differentials du^^\j = 1, A^ + 1, that is 

N 

i a (dCAdx + y^idpkA^'^P)) =0, (3.6) 

k=0 



one can determine thus the vector-functions pj = M™, j = 0, A^. As a result we 
obtain the following simple recurrent relations: 

dpj dC 



for j = 1, A^-|- 1, setting p_i = = Pn+i by definition. The unique solution to 
(3.7) is made by the following expressions, j = 0,N : 

, d'' dC 

d^dJJ+k+^) 



N ,fc 



fc=0 



Thereby we have got, owing (3.5) and (3.6), the following final representation for 
the differential dC: 

d ^ 

dC = -p[jC- J2{Pvu'^^^^^)]dx - (grad C[u],u^^^)dx + 

i=o 

d / ^ \ 

+ J- ( E (P^' ^"^'^>] + (g^^d ^M, du), (3.9) 

with -^ := di d + i d d being the Lie-derivative along the vector field 4^, and 

d^ d^ d^ ^ ^ dx' 

grad C[u\ := 6C/6u as it was mentioned above in the chapter 2. Below we intend 
to treat the representation (3.9) on the topic of a symplectic structure arisen from 
the above analysis on the invariant submanifold Mjv C M. 



4 Symplectic analysis and Hamiltonian formula- 
tion 

Let us put into the expression (3.9) the condition grad C[u\ =0 for all u = M^. 
Then the following equality is satisfied: 

, N 

dC = j-a'-^\ a(^) = J](p,-, du^^^), (4.1) 

j=o 

since the function h^^^ := '^j=q{pj, u^'-'^^^) —C{u, . . . , u^^'^^^) satisfies the condition 
dh^^^dx = —(grad C[u],u'^^^) for all x = M, owing to the relations (3.7). Taking 
now the external derivative of (4.1), we obtain that 

dx 

where we have used the well known identity d ■ -^ = -^ ■ d. From (4.2) we can 
conclude that the vector field d/dx on the submanifold Mj^ C M is Hamiltonian 
with respect to the canonical symplectic structure Vt'^'^^ = X]j=o('^Pi /\du^''^). It is a 
very simple exercise to state that the function h^^^ : J*^^^+^)(]R; M™") — > R defined 
above is playing a role of the corresponding Hamiltonian function for the vector 
field d/dx on Mjv, i.e. the equation 

dx 

holds on Mjy. Therefore, we have got the following theorem. 

Theorem 1. The critical submanifold M j^ <ZM defined by (3.1) for a given non 
degenerate smooth functional L=D(M)cD(J^^~^^^ (Il;IU")), being imbedded into 
the jet-manifold J^°°^ (Il;IU^), carries the canonical symplectic structure subject 
to which the induced vector field d/dx on M ^ is Hamiltonian. 

The theorem analogous to the that above was stated before via different man- 
ners by many authors 18, &]. Our derivation presented here is as much simpler 
as constructive, giving rise to all ingredients of symplectic theory, stemming from 
imbedding the invariant submanifold M^ into the jet-manifold. 

Now we are going to proceed further to studying the vector field (2.1) on the 
manifold Mn C M endowed with the symplectic structure f]^^) = A2( j{a^+i)(]R; ]R™))^ 
having built via the formula (4.2). 

We have the following implicating identities: 

djC dh^^^ 

^ = ^ (grad£M,irM) = -^, 

C//6 (JjJb 

dC , , ^r T du, dh^^^ , , ,. 

^ = ^ (g.ad£M,-) = -— . (4.4) 



where functions h^^' and h^^' serve as corresponding Hamiltonian ones for the 
vector fields d/dt and d/dx. This means in part that the following equations hold: 

rf/i(-) = -i^n^^\ dh^'^ = -ijL^^^l (4.5) 

dx dt 

To prove the above statement (4.5), we shall build the following quantities (for the 
vector field d/dx at first): 

d-ijL{gr ad C[u],du) = — ^(rf/i^^^) (4.6) 

d-^ dx 



stemming from (4.4), and 

■ digiad C[u], du) = -^ {ijL^^^^\ (4.7) 

dx \ d^ / 



d 



I d 



stemming from (3.9), where we used the preliminary used the evident identity 
[i_d_, ^] = 0. Adding now the expression (4.6) and (4.7) entails the following one: 

-^(grad C\ul du) = --^(dh^''^ + tjL^^^^) (4.8) 

dx dx dx 

for all a; = M and u = M. Since grad C[u] =0 for all u = M^, we obtain from (4.8) 
that the first equality in (4.5) is valid in the case of the vector field d/dx reduced 
on Mat. The analogous procedure fits also for the vector field d/dt reduced on the 
manifold Mn C M. The even difference of the procedure above stems from the 
condition on vector fields d/dt and d/dx to be commutative, [d/dt, d/dx] = what 
engenders the needed identity Ha, t"1 = hA jLi = as a simple consequence of the 

dt "^ ^dt' dx' 

procedure considered above. There upon we have stated the validity of equations 
(4.5) completely. > 

Theorem 2. Dynamical systems d/dt and d/dx reduced on the invariant sub- 
manifold Mn C M (3.1) are Hamiltonian ones with the corresponding Hamiltonian 
functions built from the equations (4.4) in the unique way. 

By the way, we have stated also that the Hamiltonian functions h^^^ and h^^^ on 
the submanifold M^ C M are commuting with each other, that is {h^^\ h^^^} = 
where {■, ■} is the Poissin structure on D{M]^) corresponding to the symplectic 
structure (4.2). This indeed follows from the equalities (4.4), since \h/~^\h'^^^} = 
^ = -^ = upon the manifold Mn C M. 

5 Symmetry invariance 

Let us consider now any vector field Kj : M ^ T{M), j = l,k, being symmetry 
fields related with the given vector field (2.1), i.e. [K,Kj] = 0, = l,k. As the 



conservation law C, G D{M) for the vector field (2.1) has not to be that for the 
vector fields Kj, j = I, k, the sub manifold M^ C M has not to be invariant also 
with respect to these vector fields. Therefore, if a vector field X G T{Mn), the 
vector field [i^j,X] ^ T{Mn) in general, if ^, j = l,k, are chosen as symmetries 

of (2.1). Let us consider the following identity for some there existing function hj : 

J(2^+2)(M; M™) ^ M, j = T/k^ stemming from the conditions [|, ^] = 0, j = T/k, 
on M : 

--i d < grad C\u\, du >= 7 — . (5.1) 

dt ^ ^ ^ dx ^ ^ 

Lemma 1. The functions hj[u], j = l,k, reduced on the invariant submanifold 
Mjy C M turn into constant ones. These constants can be chosen obviously as 
zero ones. 

Proof. We have: [^, i_d_] = 0, j = 1, k, there upon 

i _d_( i ji_d + diji) < srad C\u\,du >= -^ =^ 

dtj dt dt' dx 

=^ i_d_ I iji^d < grad C\u\, du > +di± < grad £[«], du > 

dtj \ dt dt 

dtj dt dx dtj dx dx dtj dt dx 

Whence we obtain that upon the whole jet-manifold M C J°°(]R;]R™) the 
following identities hold: 

i^(rf/iW+iAfi^^^) = /ii (5.2) 

dt j dt 

for all j = l,k. Since upon the submanifold M^v C M onr has id^Vl^'^^ = —dh'^^\ 

_ dt 

we find that hj = 0, j = 1, k, that proves the lemma. > 

Note 1. The result above could be stated also using the standard functional 
operator calculus of [8]. Indeed, 

{i^L < gradC[u],du>) (5.3) 

< grad C [u] , Kj [u] > 



dt dtj 



dt 

(^grad C[u],Kj[u])+ < grad C[u], -t^Kj[u] > 

(-K'*grad C[u],Kj[u])+ < grad C[u],K'j ■ K[u]) 
-{K'*gTa.d C[u],Kj[u])+ < grad C[u],K' ■ Kj[u]) 

d ^ ,r , ^r 1 ^^ r m dhjlu] 

-n[gmd C[u], Kj[uy - ^^ ^ 



dx dx 



Here a bilinear form 7i{-.-} is found via the usual definition of the adjoined operator 
K'* for a given operator K' : L2 -^ L2 with respect to the natural scalar bracket 



{K'*a,b):={a,K'b), {a,b) := dx{a,b), (5.4) 

Jr 

whence we simply obtain: 

{K'*a, b) - (a, K'b) = dn[a, b]/dx (5.5) 

for all a,b E L2. Therefore, we can identify hj[u\ = 7i[grad £[-u], i^j[-u]] for all 
u G M, j = l,k. li u E M]\f, we therewith obtain that hj[u] = 0, j = l,k, that 
was needed to prove. > 

As a result of the Lemma proven above one gets the following: the functions 
hj[u],j = l,k, can not serve as nontrivial Hamiltonian ones for the dynamical 
systems d/dtj,j = l,k, upon the submanifold Mn C M. To overcome this diffi- 
culty we assume the invariant submanifold M^v C M to possess some additional 
symmetries d/dtj,j = l,k, which satisfy the following characteristic criterion: 

L d grad C[u] = 0,j = 1, k, for all u G Mn- This means that for j = l,k 

dtj 

L_d_ grad C[u] = G',(grad >C[m]), (5.6) 

dtj 

where Gj{-),j = l,k, are some linear vector- valued functionals on T*{M). Other- 
wise, equations (5.6) are equivalent to the following: 

i_d_ < grad C[u], du) = —dhj[u]/dx + (7,(grad /^[-u]), (5.7) 

dt^ 



where gj{-),j = 1, A;,are some scalar linear functionals on T*{M). From (3.9) and 
(5.7) we therewith find that for all j = l,k 

L_±_ < grad C[u], du) = —-^idhj[u] + i_d_fi''^0 + dgj{gTad jC[u]). (5.8) 

If we put now u G Mjy, that is grad C[u] = 0, we immediately will find the 

following: for all j = l,k, 

dhAu]+i^n^^^ = 0. (5.9) 

dtj 

Whence we can make a conclusion that the vector fields d/dtj,j = l,k, are Hamil- 
tonian too on the submanifold Mj^ C M. Since dhj/dx = {h^^\ hj} = 0, j = 1, k, 
on the manifold M^, we therewith obtain that dh^^^dtj = 0,j = l,k. This is 
also an obvious corollary of the commutativity [d/dtj, d/dx] =0, j = 1, k, for all 
X, tj G M on the whole manifold M. Indeed, in general case we have the identity 



{h^^\ hj} = i,_d_ ^1^2*^^-*, whence the equahties {h^^\ hj} = hold for all j = 1, /c 
on the submanifold Mjv C M, since [^, ^] = on Mjy due to (5.8). The anal- 
ysis fulfilled above makes it possible to treat given vector fields d/dtj,j = l,k, 
satisfying either conditions (5.6) or conditions (5.7) on the canonical symplectic 
jet-submanifold Af/v C M analytically as Hamiltonian systems. 

6 Liouville integrability 



Now we suppose that the vector field d/dtj,j = 1, /c, are independent and commu- 
tative both to each other on the jet-submanifold Mjy C M and with vector fields 
d/dt and d/dx on the manifold M. Besides, the submanifold Mjv C M is assumed 
to be compact and smoothly imbedded into the jet-manifolcL J.*^°°^(JR; M™"). If the 
dimension dimMjsi = 2k + A, due to the Liouville theorem l[«, (jj the dynamical 
systems d/dx and d/dt are Hamiltonian and integrable by quadratures on the sub- 
manifold Mjv C M. This is the case for all Tax-mteffrable nonlinear dynamical 
systems of the Korteweg-de Vries type l[3, V, 9, 8] oil spatially one-dimensonal 
functional manifolds. 

7 Discrete dynamical systems 

Let us be given a differential discrete smooth dynamical system 

dun/dt = Kn[u] (7.1) 

with respect to a continuous evolution parameter t G M on the infinite-dimensional 
discrete manifold M C L2(Z;R'") infinite vector-sequences under the condition of 
rapid decrease in n G Z : sup„g2 1^1 ||'i*n||R'" < oo for all k G Z+ at each point 
« = (..., Un, Un+i, . . .) E M, where «„ G M™, n G Z. 

Assume further that the dynamical system (7.1) possesses a conservation law 
C G D{M), that is dC/dt = along orbits of (7.1). Via the standard oper- 
atorial analysis one gets from (3.5) the variational derivative of a functional 
^ ■= T.nez ^n[u] : 

grad Cn ■■= 6C[u]/6un = jC[*[u] ■ 1, (7.2) 

where the last right-hand operation of multiplying by unity is to be fulfilled by 
component. 

Lemma 2. Let A(M) be the infinite- dimensional Grassmannian algebra on the 
manifold M; then the differential dLn[u] G A^(M) enjoys the following reduced 
representation: 

dCn[u] =< grad £„, dun) + d/dnal^'lu], (7.3) 



where a one-form ah. [u] E A^(M) is determined in unique way, (■, ■) is the usual 
scalar product in R"^ and d/dn = A — 1, A is the usual shift operator. 

Proof. By definition we obtain for tlie external differential dCn[u\ the following 
chain of representations for each n G Z : 



N 



dCn[u] =^< 



dCn[u\ 



k=0 



du. 



n+k 



dUn+k > 



>r^>r;^ (i dCn~s[u] , , >r^ dCn-k[u] , 



dn dun+k-. , „ 

fc=0 s=0 fc=0 



d 

dn 



N k 



dCn-Ju] 



N 



X] X^ < :i^-^^^i^, dUn+k-s > 



k=0 s=0 



du 



n+k—s 



k=-N 



< 



dCJu] 



du 



,dUn > 



n+k 



d v^ v^ dCn-s[u] ^..^ 

-r2^2^< -^. ' dUn+k-s > + < C^ .l,dUn> 



dn 



k=0 s=0 



du 



n+k—s 



d_ 

dn 



al>[u]+ < grad Cn[u],dUn), 



(7.4) 



where N G Z^ is the fixed number depending on the jet-form of a functional 

CeD{M), 

dCn-s [u] 



al'^Ml 



N k 

EE 



< 



fc=0 s=0 

TV k 



du 



n+k—s 



-, dUn+k-s > 



EE<^%^. ''""«>. 



A;=0 jr=0 



9-u^ 



(7.5) 



n+j 



and 



jrad £„[m] = £„ -1 = 2^ ^ • 

fe=0 " 



The latter equality in (7.4) proves the lemma 2 completely. > 

The above proved representation (7.3) gives rise to the following stationary 
problem being posed on the manifold M : 



Mn = {u e M : grad £„ = 0} 



(7.6) 



for all n G Z, where by definition det 
the validity of the following theorem. 



d'^Cnlu] 



9«]V+i(9mjv+i 



0. In virtue of (7.3) we obtain 



Theorem 3. The finite- dimensional Lagrangian suhmanifold Mn C M defined 
by (7.6), is a symplectic one with the canonical symplectic structure r2„ = dan 
being independent of the discrete variable n G Z. 

Proof. From (7.3) we have that on the manifold M^ C M dCn[u] = d/dn{an [u]), 
whence for all n E Z d/dn{fln ) = 0. This means obviously, that fi„+i = fin 
for all n G Z, or equivalently, the 2-form fin is not depending on the discrete 
variable n G Z. As the 2-form fi„ := dali by definition, this form is chosen to be 
a symplectic form on the manifold Mn C M. For this 2-form to be nondegenerate 



d^Cnlu] 



dUn+N+ldUn+N+1 



^ on M, 



TV- 



on Mjv, we assume that Hessian of £„ equals det 
The latter proves the theorem. [> 

Let us consider now the prior given dynamical system (7.1) reduced on the 
submanifold Mj^ C M. To present it as the vector field d/dt on M^v, we need 
preliminary to represent it as a Hamiltonian flow on M^y. To do this, let us write 
the following identities on M : 

id_d < grad Cn,dun) = — —iA^ni^], 
dt dfi dt 

di± < grad £„, (iMn) = — 3— (c^^n M), (7.7) 

•^t dn 

which are valid for all n G Z. Adding the last identities in (7.7), we come to the 
following one for all n G Z : 

I < grad £„,rfw„) = -A(^^fi(f)M +/^W[^]). (7.8) 

Having reduced the identity (7.8) upon the manifold Mjv C M, we obtain the 
wanted expression for all u G Mat, N E Z : 

^.fi(f)M + /^WM = 0. (7.9) 

dt 

The latter means that the dynamical system (7.1) on the manifold Af/v is a Hamil- 
tonian one, with the function /in [m] being a Hamiltonian one defined explicitly by 
the second identity in (7.7). 

We assume now that the symplectic structure fin [u] on M^v be representable 
as follows: 

N 

^f N = 5^ (rfPi+n A dUj+n) , (7. 10) 

i=0 



where generalized coordinates Pj+n ^ IR™, j = 0, N, are determined from the 



following discrete jet-expression Cn{u\ := £(«„, Wn+i, • • • , Un+N+i)^ n E7i, 

i=o fc=o j=o "'+■5' 

^ ^ (9/: _ 

3=0 k=j "+■' 



whence we get the final expression: 



dCn+j-k[u] 






fc=i 



where j = 0, iV, m G Mat C M. 

Now we are in a position to reformulate the given dynamical system (7.1) as 
that on the reduced manifold M^ C M : 



for all n E Z,j = 0,N. Thereby the problem of embedding a given discrete dy- 
namical system (7.1) into the vector field fiow on the manifold Mjv C M is solved 
completely with the final result (7.12). 

8 Invariant Lagrangians construction: functional 
manifold case 

In the case the given nonlinear dynamical system (2.1) being integrable one of 
Lax-type, we can proceed effectively to find a commuting infinite hierarchy of 
conservation laws can be serving as the invariant Lagrangians, having been under 
consideration above. 

At first we have to use the important property l[3] of the complexified gradient 
functional ip = grad 7 G T*{M) ® C generated by an arbitrary conservation law 
7 G D{M), i.e. the following Lax-type equation : 

dip/dt + K'*ip = (8.1) 

where the prime sign denotes the usual Frechet derivative of the local functional K : 
M — > T{M) on the manifold M, the star "*" denotes its conjugation operator with 
respect to the nondegenerate standard convolution functional (■, ■) = f^dx{-, ■) on 



T*{M) X T{M). The equation (8.1) admits (what follows from |3U, 3i, 32]) the 
special asymptotic solution: 

(fix, t; A) = (1, a{x, t; X)yexp[uj{x, t; A) + d~^a{x, t; A)], (8.2) 

where a{x,t;\) G C™""*^, a{x,t;\) G C, uj{x,t;X) - some dispersive function. 
The sign "r" means here the transposition one, what is adopted in matrix analysis. 
For any complex parameter A G C at |A| ^ oo the following expansions take place: 

a{x,t;X) ~ J2 ciA^,t;u]X-^+''^^\ a{x,t;X) ~ ^ aj[x,t;u]X-^+''^''\ 
iez+ jez+ 

Here s{a) and s{a) G Z_|_ - some appropriate nonnegative integers, the operation 
d^^ means the inverse one to the differentiation d/dx, that is d/dx ■ d^^ = 1 for 
all X eR. 

To find the explicit form of the representation (8.2) in the case when the associ- 
ated Lax-type representation |»] depends parametrically on the spectral parameter 
A(t; z) G C, satisfying the following non-isospectral condition: 

d\{t;z)/dt = g{t;\{t;z)), A(t; 2)|^^o+ = ^ G C, (8.3) 

for some meromorphic function g{t; ■) : C ^ C, t G IR+, we must reanalyze more 
carefully the asymptotic solutions to the Lax equation (8.1). Namely, we are going 
to treat more exactly the case when the solution (f G T*{M) to (8.1) is represented 
as an appropriate trace-functional of a Lax spectral problem at the moment r = 
t G M+ with the spectral parameter A(t; A) G C satisfying the condition (8.3), the 
evolution of the given dynamical system (2.1) being considered with respect to the 
introduced above parameter r G M, that is 

du/dr = K[x, r; u], (8.4) 

u\^^Q = u & M - some Cauchy data on M. This means that the functional 

ip{x,T;\) := reggia.d Sp S{x,t;X), \ = \{t; \(t; z)) E C, (8.5) 

where S{x, r; A) is the monodrony matrix corresponding to a Lax type spectral 
problem assumed to exist, has to satisfy the corresponding Lax equation at any 
point u E M subject to (8.4): 

dl^/dT + K'*[u]-l^ = Q (8.6) 

for all r G M+. Under the above assumption it is obvious that the spectral param- 
eter A = A(r; A(t; z))^ where 

d\/dT = g(T;~X), A =\(t;z)EC, (8.7) 

r=0 



g{t; ■) : C ^ C - some meromorphic function being found simply from (8.6) for 
instance at m = 0, the Cauchy data A(t; z) G C for all t G M+ are the corresponding 
ones to (8.3), the parameter z G C being a spectrum value of the associate Lax 
type spectral problem at a moment t G M+. 

Now we are in a position to formulate the following lemma. 

Lemma 3. The Lax equation (8.6) as the parameter r = t E M+ admits an 
asymptotic solution in the form 

^(x, r; A) = (1, a{x, r; \)yexp[uj{x, r; A) + d^^a^x, r; A)], (8.8) 

where d{x, r; A) G C"'''^, o-{x, r; A) G C, are some local Junctionals on M, u{x, r; A) G 
C is some dispersion function for all x eM,t E IR+, and if for |A| ^ oo the prop- 
erty \\\ ^ oo as T = t E IR+ holds, the following expansions follows: 

a{x,r;\) ~ >^ aj[x,T;u]\ , a(x, r;A) ~ >^ ^^[a;, r; -ujA , (8.9) 

iez+ jez+ 

with s{a) and s(a) G Z+ being some integers. 

Proof. In virtue of the theory of asymptotic expansions for arbitrary differential 
spectral problems, the result (8.8) will hold provided the representation (8.5) is 
valid and the spectral parameter \(t; z) G C is taken subject to (8.7). But the 
above is the case because of the Lax-type integrability of the dynamical system 
(8.4). Further, due to the mentioned above integrability of (8.4) as well as due 
to the well known Stokes property of asymptotic solutions to linear equations like 
(8.1), the condition (8.3) holds for some meromorphic function g{t; ■) : C — > C, 
t G M+, enjoing the determining property -^ J^a{x,t; \{t; \(t; z)))dx = for all 
t G M+. The latter proves the lemma completely. > 

As a result of Lemma 3 one ca n formulate the following important theorem. 

Theorem 4. The Lax integrable parametrically isospectral dynamical system 
(8.4) as T = t E M+ admits an infinite hierarchy of conservation laws, in gen- 
eral nonuniform ones with respect to the variables x E M.,t E M+, which can be 
represented in an exact form in virtue of the asymptotic expansion (8.8) and (8.9). 

Proof. Indeed, due to the expansion (8.8), we can obtain right away that the 
functional 

^{t]\(t]z))= dxa{x,t;\{t;\{t;z))) (8.10) 



does not depend at r = t G M+ on the parameter t E M_|_, that is 

d^/drl^,^^^ = (8.11) 

for all t E M+. If we put also the parameter r G M+ to tend to t G IR+, due 

to (8.5) we obtain that ip{x,T;X) — * ip{x,t;X) for all a; G M, t G M+ and 

T=teR+ 



A(t; z) G C This means that a complexified local functional ip{x, t; z) G T*(M) ®C 
satisfies the equation (8.1) at each point u G M. As an obvious result, the following 
identifications hold: 



a;(x, r; A) 



=teK+ 



uj{x,t;z), a{x,T;X) 



a(x,t: z) 



=teiR4 



for all 2; G C. Hence, the functional 

7(z) := 7(r; A(t; 2;))|^^jg]g = Jj^ (ia;(T(a;, t; z) G D{M) doesn't depend on the evo- 
lution parameter t G M+ and due to equation (8.1) is a conserved quantity for the 
nonlinear dynamical system (2.1) under consideration, i.e. 



d-f{t;z)/dt = 



^.12) 



for all t G IR+ and 2;' G C. Therefore, it makes it possible to use the equation 
(8.12) jointly with (8.7) for the asymptotic expansions (8.9) and (8.3) to be found 
in exact form. To do this we at first need to insert the asymptotic expansion (8.8) 
in the determining equation (8.6) for the asymptotic expansions (8.9) to be found 
explicitly at the moment t = t E M+. Keeping in mind that at r = t G M+ |A| -^ 
00 if I A| -^ 00, and solving step by step the resulting recurrence relationships for the 
coefficients in (8.9), we will get the functional 7(2;) := 7(r; A(t; z))\^^^^-g^ , z G C, 
in the form fitting for the criteria equation (8.12) could be used. As the second 
step, we need to use the differential equation (8.7) for the criteria equation (8.12) 
to be satisfied point- wise for all t G M+. This means, in particular, that 



dl{z) 
dt 






dt 



ieZ4 



dx(yj\x^ r; -ulA 



~ -i+s(5-) 



(8.13) 



dx y 
dx y 

iez+ 



~di +ctj[x,t;u]X (^(^)~-^)^ 



T=t£R+ 



{daj/dt)X 



~-j+s{o-) 



A;>>— 00 



{s{a) - k)akgj-k-i{t)\ 



+ 



iGZ+ 



i.s{(y)~3)^g{t;\) 



r-fe 



where we have put by definition g{T;\) :~ '}2k»~oo9k{''~)-^ fo^ ''" ^ ^1^+ ^^^ 
|A| ^ 00. Since the spectral parameter A = A(t; z) at the moment t = 0^ coincides 



with an arbitrary complex value 2; G C, the condition |2;| ^ 00 together with (8.13) 
at the moment t = 0^ gives rise to the following recurrent relationships: 



E 

iGZ+ 



daj/dt + a^- K[t;u]+ ^ 



s o- 



A;)o-fc • Qj^k- 



k»-c 



~-j+s(d-) 



=tgR+ 



(8.14) 



dX 



Yl ^ji-^i^)-j)j^9it;\)\ 



~-j+s(a)-i 



iGZ4 



{mod d/dx) 



=teR+ 



for all j G Z+, x EM.^t E M+ and u E M. Having solved the algebraic relationships 
(8.14) for the prior unknown function g{t; X),t G M+, we will obtain the generating 
functional 7(2;), 2; G C, of conservation laws for (2.1) in exact form. This completes 
the constructive part of the proof of the theorem above. 

From the practical point of view we need first to get the differential equation 
(8.7) in exact, maybe in asymptotic form and find further the dispersive function 
ij{x,t; A) and the local generating functional d'{x,T; A) defined via (8.8) and (8.9) 
for all a; G M, r G M+ and |A| -^ 00, and next one can find the equation (8.3) due 
to the algorithm based on the relationships (8.14). This together with the pos- 
sibility of applying the general scheme of the gradient-holonomic algorithm 1 50] 
gives rise to determining in many cases the above mentioned Lax-type represen- 
tation completely in exact form, that successfully solves the pretty complex direct 
problem of the integrability theory of nonlinear dynamical systems on functional 
manifolds. 

Having obtained the generation function '-/{z) G D{M), z E C, oi an infinite 
hierarchy of conservation laws of the dynamical system (2.1) on the manifold M, 
we can build appropriately a general Lagrangian functional Cn ^ D{M) as follows: 



N 



c 



N 



-7Jv^ 



where, by definition, 

^^dxaj[x,t]z\, j G Z+, 



j=0 

7(z) = f^dxa{x,t; z) and for 



(8.15) 



00 functionals 7- 



are conservation laws due to expansion (8.2), with Cj G 



0, A^, being some arbitrary constants and A^ G Z+ being an arbitrary 
but fixed nonnegative integer. If the differential order of the functional 7jv+i ^ 
D[M) has the highest one of the orders of functionals 7- G D{M),j = 0,N, and 
additionally, this Lagrangian is not degenerate, that is det{Hess'y^^i) 7^ 0, we can 
apply in general amost all the theory developed before, to prove that the critical 
submanifold M^ = {u E M : grad Cn = 0} is a finite-dimensional invariant 
manifold inserted into the standard jet-manifold J*^°°^(M;M™) with the canonical 
symplectic structure subject to which our dynamical system is a finite-dimensional 
Hamiltonian flow on the invariant submanifold M^. 



9 Invariant Lagrangian construction : discrete 
manifold case 

Let us consider the discrete Lax integrable dynamical system on the discrete mani- 
fold M without an a priory given Lax-type representation. The problem arises how 
to .get the corresponding conservation laws via the gradient-holonomic algorithm 
ls]7 To realize this way let us study solutions to the Lax equation: 

d^Jdt + K[r,u]-^^ = Q (9.1) 

local functional v^„[m] G T*^{M) at the point «„ G M, n G Z. In analogy to 
the approach of Chapter 7 we assert that equation (9.1) admits a comlexified 
generating solution (p^ = (pj^it; X) G T*^{M) ® C, n G Z, with z E C being a 
complex parameter in the following form: 

<f,{t- z) ^ (1, an{t; z)yexp[u{t- z)] ( J] a,{t; z)\ , (9.2) 

\;=— oo ^ 

where cj(t; z) is some dispersive function for t G M+, a„(t; z) G C'"^"'^, 

o"„(t; 2;) G M are local functionals on M, having the following asymptotic expansions 

at \z\ -^ 00: 

a„(t; z)c^Yl «"[^' ^]^''^-'^''\ ^i'H^; ^)-Yl ^"[^' u]^''^'"^"^- (9-3) 

iGZ+ jGZ+ 

To find the explicit form of the asymptotic representation (9.2) we need to study 
additionally the asymptotic solutions to the following attached Lax equation with 
respect to the new evolution parameter r G M+: 

dipJdT + K'^*[T,u]-^^ = 0, (9.4) 

where (^„ G T*^{M) ® C, and a point u E M evolves subject to the following 
dynamical system: 

dun/dr = Kn[T; u], (9.5) 

for all n G Z. Having made the assumption above we can assert based on the 
general theory of asymptotic solutions to linear equations like (9.4), that it admits 
also in general another asymptotic solution in the similar form: 

n 

^„(r; A) ^ (1, o„(r; A))^exp[cu(r; A)] J] ^.(^S ^)' (9-6) 

i=-oo 



where for all n G Z and at r G IR+ the asymptotic expansions 

iez+ 

i€Z+ 



(9.7) 



hold. The expansions above are valid if |A| — ;► oo as |A(t;z)| — > cxd, z g C The 
latter is the case because of the Lax-integrability of the dynamical system (9.5). 
The evolution 



dX/dr = gij] A), A 



t=o 



\{t]z) gC, 



(9.^ 



where ^(r; •) : C ^ C is some meromorphic mapping for all r G M+, is in general 
found making use of the corresponding solution to (9.4) at -u = 0. 

Substituting the expansions (9.6) and (9.7) into (9.4), we obtain some recur- 
rence relationships, giving rise to a possibility the expressions for local functionals 
5"j[t;M„],j G Z+, which can be found exactly. Having this done successfully, we 
assert that the functional 



lit'^z) = ^ In o-„(r; A) 



^lno-„(t; 



(9.9) 



neZ neZ 

where A = A(r; A), r G IR+, and A(t; z) G C, is a meromorphic solution to equation 



d\/dt = g{t;\), A| 



t=o+ 



zeC 



(9.10) 



with still independent meromorphic function g{t; ■) for almost all t G M+. The 
latter can be found making use of the following determining condition: the local 



functional V'„(t; A) 



(Pn{t;z) G T*{M) ® C for all t G M+ and z G C. 



Hence, the following equality holds immediately: 

— (^ lnCT„(r;A) 



dt 



E ^n\t;X) 



nGZ4 



nGZ+ 

da 



(9.11) 



da.. 



^^ +a'^.K4u] + -^git;\) + 



dan d\ , , , 
dX dX""^ 







for all t G IR+. Equating coefficients of (9.11) at all powers of the spectral parameter 
A(t; 2;) G C to zero modulus d/dn, n G Z, we will find the recurrent relationships 



for the function g(t; A) of (9.8) to be determined successfully. Thereby, using the 
equation (9.10) and an expansion a{t;z) ~ 12jez '^j[^jUn\z~^~^^^^^ for \z\ -^ cxd, 
where s{a) G Z+ is some integer number, we obtain an infinite hierarchy of discrete- 
wise conservation laws of the initially given nonlinear dynamical system (2.1) on 
the manifold M. But because of the parametric dependence of the conservation 
laws built above on the evolution parameter t G M+, we cannot use right now the 
theory developed before to prove the Hamiltonian properties of the corresponding 
vector fields on the invariant submanifolds. To do this in an appropriate way, it is 
necessary to augment the theory developed before in some important details. 

10 The reduction procedure on nonlocal Lagrangian 
submanifolds 

1. The general algebraic scheme. Let Q := C°°(E>^;Q) be a Lie algebra of 
loops, taking values in a matrix Lie algebra Q. By means of Q one constructs the 
Lie algebra Q of matrix integral-differential operators \'3^: 

a:=Y, «.e> (10.1) 

j<Coo 

where the symbol ^ := d/dx signs the differentiation with respect to the inde- 
pendent variable x G M/27rZ ~ §^. The usual Lie commutator on Q is defined 

as: 

[d,b] := aob-bod (10.2) 

for all d,b & Q, where "o" is the product of integral-differential operators, taking 
the form: 

aez+ ^ 

On the Lie algebra Q there exists the adinvariant nondegenerate symmetric bilinear 
form: 

(a, 6) := / Tr{dob) dx, (10.4) 

Jo 

where Tr-operation for all a G ^ is given by the expression: 

Trd:=res^Spd = Spa^i, (10.5) 

with Sp being the usual matrix trace. With the scalar product (10.4) the Lie 
algebra Q is transformed into a metrizable one. As a consequence, the dual to Q 



linear space of the matrix integral-differential operators Q* is naturally identified 
with the Lie algebra Q, that is Q* ~ Q. 

The linear subspaces Q+ G Q and Q^ G Q such as 

n(a)<?;oo 

g+:= <j a := ^ ttj^^ : ttj eGJ = 
i=o 



oo 




^-^= 1 ^ •= E^'^'^'^^^- ^ ^^ e ^, J G Z+ ^ , (fO.6) 

I i=o 

are Lie subalgebras in Q and Q = Q+ ® Q-- Because of the splitting of Q into 
the direct surn of its Li,ie subalgebras one can construct the so called Lie-Poisson 
structure p/, 5«', 5y,' 4u] on Q* , using a special linear endomorphism 7^ of ^: 

7^:=(P+-P_)/2, P±^:=^±, P±^^ = 0. (10.7) 

For any smooth by Frechet functionals 7,/i G V{Q*) the Lie-Poisson bracket 
on Q* is given by the expression: 

{7, ^}^ (0 = (/, [V7(/), Vmn) , (10.8) 

here^ / G Q* and for all a, b G Q the 7?.commutator in (10.8) has the form 
^750]: 

[a, % := [TZa, b] + [a, TZb], (10.9) 

subject to which the linear space Q becomes a Lie algebra too. The gradient 
V7(/) G ^ of a functional 7 G Vi^Q*) at a point / G ^* with respect to the scalar 
product (10.4) is defined as 

<57(/):=(V7(0,5/^), (10.10) 

where the linear space isomorphism ^ ~ ^* is taken into account. 

The Lie-Poisson bracket (10.8) generates Hamiltonian dynamical systems on 
Q* related with Casimir invariants 7 G I{Q*), satisfying the condition: 

[V7(/),/] = 0, (10.11) 

as the corresponding Hamiltonian functions. Due to the expressions (10.8) and 
(lO.ll)the mentioned above Hamiltonian system takes the form: 

dl/dt := [7^V7(0, i] = [V7+(0, ^1, (10.12) 



I PM La 

being equivalent to the usual commutator Lax type representation l[4u, 4i]. The 
relationship (10.12) is a compatibility condition for the linear integral-differential 
equations: 

if = A/, 

df/dt = V7+(/)/, (10.13) 

where A G C is a spectral parameter and a vector-function / G W^(S"'^;H) is an 
element of some matrix representation for the Lie algebra Q in some functional 
Banach space H. 

Algebraic properties of the equation (10.12) together with (10.14) and the 
associated dynamical system on the space of adjoint functions /* G H^*(S^; H): 

rfrM = -(V7(/));r, (io.i4) 

where /* G W* is a solution to the adjoint spectral problem: 

i*f* = uf*, (10.15) 

being considered as some coupled evolution equations on the space Q* (BW Q) W* 
is an object of our further investigation. 

2. The tensor product of Poisson structures and its Backlund trans- 
formation. To compactify the description below we will use the following desig- 
nation of the gradient vector 

V7(/;/,f) := {6i/6l, 6^/6f, 6^/6rf 

for any smooth functional 7 G V{Q* © WfB tF*) On the spaces Q* and W © W* 
there exist canonical Poisson structures |3», 42, 40] 

S^/Sl ;i. [{5^/5lU, I] - [6^/51, lU (10.16) 

at a point / G Q* and 

{6^/6 f, 6^/6rf -A {6^/6r, -6^/6 ff (10.17) 

at a point (/, /*) E W (B W* correspondingly. It should be noted that the Poisson 
structure (10.17) is transformed into (10.12) for any Casimir functional 7 G I{Q*) 
. Thus, on the extended space ^* © W^ © W* one can obtain a Poisson structure 
as the tensor product 6 := 6^ © J of the structures (10.17) and iJO.lS). 
Let us consider the following Backlund transformation I4u, 42, 4b]: 

{I, f, n -A {i{i, f, n,f = f, r = n, (lo.is) 



generating on Q* Q)W Q)W* a Poisson structure 9 with respect to variables (/, /, /*) 
of the coupled evolution equations (10.12), (10.14), (10.15). 

The main condition for the mapping (10.19) to be defined is the coincidence of 
the dynamical system 

{di/dt, df/dt, df/dtf := -eV7(/, /, /*) (10.19) 

with (10.12), (10.14), (10.15) in the case of 7 G I{G*), i.e. if this functional is taken 
to be not dependent of variables (/, /*) E W Q) W*. To satisfy that condition, one 
has to find a variation of any smooth Casimir functional 7 G I{Q*) as 61 = 0, 
considered as a functional on Q* (B W (B W*, taking into account flows (10.14), 
(10.15) and the Backlund transformation (10.19): 

siil f, h ,- = (< h/sL sf >) + (< s^/sr, 6r >) = 
(< -dr/dt,5f >) + (< df/dt,5r >)\.^^ ^-^^^^ = 

(< {6^/6irj*,6f >) + (< {6^/6l)^f,6r >) = 
(< /*, (57/50+5/ >) + (< (57/5/)+/, 5/* >) = 

(57/5/, 6fr' ® /*) + (57/5/, fc' ® 5/*) = 

(57/5/, 5(/r' ® /*)) := (57/5/, 51). (10.20) 

As a result of the expression (10.21) one obtains the relationships: 



61 



= 5(/r'®r), (10.21) 



or having assumed the hnear dependence of / and / G Q* one gets right away that 

i=i + fr'®r- (10.22) 

Thus, the Backlund transformation (10.19) can be now written as 

{ijj*)-A{i=i-fc'®r, f = f, r = h- (10.23) 

IPM SP 

The expression (10.24) generalizes the result, obtained in the papers |4Li, 46] for 
the Lie algebra Q of integral-differential operators with scalar coefficients. The 
existence of the Backlund transformation (10.19) makes it possible to formulate 
the following theorem. 

Theorem 1. A dynamical system on Q* (B W (B W* , being Hamiltonian with 
respect to the canonical Poisson structure : T* {Q* Q)W Q)W*) -^ T{g*®W®W*), 
and generated by the evolution equations: 

dl/dt = [V7+(0, /] - [V7(0, lU, df/dt = 5^/5r, dr/dt = -6^/6 f, (10.24) 



with 7 G I{Q*) being the Casimir functional at I E Q* connected with I E Q* by 
(10.23), is equivalent to the system (10.12), (10.14) and (10.15) via the constructed 
above Backlund transformation (10.24). 

By means of simple calculations via the formula (see for egz. l[4Li, '66]) 

e = b'qb'*, 

where B' : T{g* ®W® W*) -^ T{g* ®W® W*) is the Frechet derivative of (24), 
one brings about the following form of the Poisson structure O on ^* © W^ © W* 3 

(/,/,/*): 

\-fr' ® ^7/5/ + s^/sf*r' ® /* sj/sf* - (<57/5/)+/ - sj/6f + {sj/si)xf, 

that makes it possible to formulate the next theorem. 

Theorem 2. The dynamical system (10.20), being Hamiltonian with respect to 
the Poisson structure in the form (10.26) and a functional 7 G I{G*), gives the 
inherited Hamiltonian representation for the coupled evolution equations (10.12), 
(10.14), (10.15). 

By means of the expression (10.23) one can construct Hamiltonian evolution 
equations, describing commutative flows on the extended space ^* © W^ © W* at 
a fixed element I E Q* . Due to (10.24) every equation of such a type is equivalent 
to the system 

di/dTn=[llX 

df/dTr, = ilf, (10.25) 

df*/dT„ = -ihir, 

generated by involutive with respect to the Poisson bracket (10.17) Casimir invari- 
ants 7„ G I{g*), n G N, taking here the standard form: 

7„ = l/(n+l)(r,/) 

at / G Q*. 

The compatibility conditions of the Hamiltonian systems (10.25) for different 
n G Z_(_ can be used for obtaining Lax integrable equations on usual spaces of 
smooth 27r-periodic multivariable functions that will be done in the next section. 

3. The Lcix type integrable Davey-Stewartson equation and its triple 
linear representation. Choose the element / G ^* in an exact form such as 

' I -1 r [ u 



where u,u e C~(§^ C) and G = gl{2; C). Then 

where / = {fi, f-z)^ and /* = (/i*,/2)"^, "~" can sign the complex conjuga- 
tion. Below we will study the evolutions (10.25) of vector-functions (/, /*) G 
iy(§i;C2) © W^*(§i;C2) with respect to the variables y = n and t = ra at the 
point (10.26). They can be obtained from the second and third equations in 
(10.25), having put n = 1 and n = 2, as well as from the first one. The latter is 
the compatibility condition of the spectral problem 

/$ = A$, (10.27) 

where $ = ($i,$2)"^ £ W^(§"'^;C^), A G C is some parameter, with the following 
linear equations: 

d^/dy = /+$, (10.28) 

d^/dt = /^$, (10.29) 

arising from (10.26) at n = 1 and n = 2 correspondingly. The compatibility of 
equations (10.28) and (10.29) leads to the relationships: 

du/dy = -2/i/*, du/dy = -2flU, (10.30) 

dfi/dy = dfi/dx - w/a, dfl/dy = dfl/dx - uf;, 

df2/dy = -dh/dx + M/i, df;/dy = -df;/dx + ufi 

Analogously, replacing t G C by it G iM, i^ = —1, one gets from (10.29) and 
(10.30): 

du/dt = t{d\/dxdy + 2u{hfl + /2/2)), du/dt = -i{d^u/dxdy + 2u{hfl + /2/2)), 
dUifD/dy - dU\n)/dx = l/2d{uu)/dy = -{d{M;)/dx + d{hf;)/dy), 
dh/dt = t{d'f\/dx' + {2j\f* - uu)h - du/dxh), 

df*/dt = -tid'rjdx' + (2/i/* - uu)r, - du/dxf;), (10.31) 

df2/dt = i{d^h/dx^ - (2/2/2* + nu)f2 - du/dxfi), 

df;/dt = -i{d^f;/dx^ - (2/2/* + uu)f; - du/dxf*). 



The relat 
equation 



(iishim (10.31), (10.32) take the well known form of the Davey-Stewartson 



3», 43] at u gC°°(§^;C) being a complex conjugated to u gC°°(S 



oo/sl. 



The compatibility for every pair of equations (10.28), (10.29) and (10.30), which 



can be rewritten as the first order linear ordinary differential ones in such a way: 

A u -/i 
d<l>/dx = \ u -X /2 I $, (10.32) 



d<^/dy =0 A -/2 $, (10.33) 



U -A /2 

/r /2 ^ 


j*' 






A -/i \ 
A -/2 

/r /2 / 


<^, 






A' + /i/r 

-1/29m/9i/ 


1/29m/9i/ 

A' - /2/2* 


-A/i- 
-A/2- 


- dh/dy 



ci<l>/cit = t I -l/29iz/9i/ A^ - M; -a/2 - 9/i/9|/ I $(10.34) 

A/r + 5/i75y \f; + df;/dy o 

where $ = ($1, $2, ^^s)"^ ^ W^(§^;C^), provide its Lax type integrability. Thus, 
the following theorem holds. 

Theorem 3. The Davey- Stew arts on equation (10.32), (10.33) possesses the 
Lax representation as the compatibility condition for equations (10.34) under the 
additional natural constraint (10.27). 

In fact, one has found above a triple linearization for a (2+l)-dimensional dy- 
namical system, that is a new important ingredient of the Lie algebraic approach 
to Lax type integrable flows, based on the Backlund type transformation [10.23) 
developed in this work. It is clear that the similar construction of a triple lineariza- 
tion like (10.4) can be done for many other both old and new (2+l)-dimensional 
dynamical systems, on what we plant to stop in detail in another work under 
preparation. 

11 Conclusion 

The developed above theory of parametrically Lax-type integrable dynamical sys- 
tems concedes to widen to a great extent the class of exactly treated nonlinear 
models in many fields of science. It is to be noted here the following impor- 
tant mathematical fact being got in the paper: almost every nonlinear dynamical 
system admits a parametrically isospectral Lax type representation but a given 
dynamical system is the Lax-type integrable if an evolution of the spectrum pa- 
rameter doesn't depend on a point m G M at all Cauchy data. This result has 
allowed us to develop a very effective direct criterion for the following problem: 
whether a given nonlinear dynamical system on the functional manifold M is para- 
metrically Lax- type integrable or not. Having the problem above solved, we have 
suggested the reduction procedure for the associated nonlinear dynamical systems 
to be descended on the invariant submanifold M^v C M built before inheriting the 
canonical Hamiltonian structure and the Liouville complete integrability. Thereby, 



the powerful techniques of perturbation theory can be successfully used for dy- 
namical systems under consideration, as well as the relationships between the full 
Hamiltonian theory and various Hamiltonian truncations could be now got under- 
standable more deeply. 

The imbedding problem for infinite-dimensional dynamical systems with ad- 
ditional structures such as invariants and symmetries is as old as the Newton- 
Lagrange mechanics, having been treated by many researches, using both analyt- 
ical and algebraic methods. The powerful differential-geometric tools used here 
were created mainly in works by E. Cartan at the beginning of the twentieth cen- 
tury. The great impact in the development of imbedding methods was done in 
last time, especially owing to theory of isospectral deformations for some linear 
structures built on the special vector bundles over the spase M as the base of a 
given nonlinear dynamical system. Among them there are such structures as the 
moment map I : M —^ Q* into the adjoint space to the Lie algebra Q of symmetries, 
acting on the symplectic phase space M equivariantly |8, lu], the connection of 
the Cartan-Eresman structures appearing via the Wahlquist-Estabrook approach 
[TT], and many others. 

For the last years the general structure of Lagrangian and Hamiltonian for- 
malisms was studied thoroughly using both geometrical and algebraical methods 
I i2, 13]. The special attention was paid to the theory of differential-difference dy- 
namical systems on the infunite-dimensional manifolds 1 13, ?T. Some number of ar- . 
tides was devoted to the theory of pure discrete dynamical systems l[ib, 16, lY, i«], 
as well treating the interesting examples Ifl^] appeared to be important for appli- 
cations. 

In future work we intend to treat further imbedding problems for infinite- 
dimensional both continuous and discrete dynamical systems basing on the differential- 
geometric Cartan's theory of differential ideals in Grassmann algebras over jet- 
manifolds, intimately connected with the problem under regard. As it is well 
^npwm there existed by now only two regular enough algorithmic approaches 
\'6()' 4u,'46] to constructing integrable mult i- dimensional (mainly 2+1) dynami- 
cal systems on infinite-dimensional functional spaces. Our approactL„deyised in 
this work, is substantially based on the results previously done in l[4(j, 46], ex- 
plains completely the computational properties of multi-dimensional flows before 






delivered in works 



'6'6, 45]. As the key points of our approach there used the 
canonical Hamiltonian structures naturally existing on the extended phase space 
and the related with them Backlund transformation which saves Casimir invari- 
ants of a chosen matrix integral-differential Lie algebra. The latter gives rise to 
some additional Jlaroiltonian properties of considered extended evolution flows be- 
fore ^udied in l[5«, 40] making use of the standard inverse scattering Jg-ai^orm 
141', 5o, 40] and the formal symmetry reduction for the KP-ierarchy \Ai^ 46] of 



commuting operator flows. 

As one can convince ourselves analyzing the structure of the Backlund type 
transformation (10.24), that it strongly depends on the type of an arf-invariant 
scalar product chosen on an operator Lie algebra Q and its Lie algebra decompo- 
sition like (10.6). Since there exist in general other possibilities of choosing such 
decompositions and ad-invariant scalar products on Q^ they give rise naturally to 
another resulting types of the corresponding Backlund transformations, which can 
be a subject of another special investigation. Let us here only mention the choice 
of a scalar product related with the operator Lie algebra Q centrally extended by 
means of the standard Maurer-Cartan two-cocycle l|36,'3y, '40j, bringing about new 
types of multi-dimensional integrable flows. 

The last aspect of the Backlund approach to constructing Lax type integrable 
flows and their partial solutions whichis worth of mention is related with Darboux- 
Backlund type tainsformations p^, 44] and their new generalization recently de- 
veloped in |4b', 46]. They give rise to very effective procedures of constructing 
multi-dimensional integrable flows on functional spaces with arbitrary number of 
independent variables simultaniously delivering a wide class of their exact analyt- 
ical solutions, depending on many constant parameters, which can appear to be 
useful for diverse applications in applied sciences. 

All mentioned above Backlund type transformations aspects can be studied 
as special investigations, giving rise to new directions in the theory of multi- 
dimensional evolution flows and their integrability. 
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